A VARIATIONAL APPROACH TO THE MACROSCOPIC 
ELECTRODYNAMICS OF ANISOTROPIC HARD 
SUPERCONDUCTORS 



GRAZIANO CRASTA AND ANNALISA MALUSA 

Abstract. We consider the Bean's critical state model for anisotropic superconductors. 
A variational problem solved by the quasi-static evolution of the internal magnetic field 
is obtained as the F-limit of functionals arising from the Maxwell's equations combined 
with a power law for the dissipation. Moreover, the quasi-static approximation of the 
internal electric field is recovered, using a first order necessary condition. 

If the sample is a long cylinder subjected to an axial uniform external field, the 
macroscopic electrodynamics is explicitly determined. 



1. Introduction 

It is well known that a superconductor is a conductor which is able to pass an elec- 
tric current without dissipation. The transition from the normally conducting state to 
the superconducting one occurs at a critical temperature T c , depending on the material, 
below which the material exhibits (almost) perfect conductivity. We are interested in the 
response of a superconducting material to an applied external magnetic field H s under 
isothermal conditions below its critical temperature. 

The superconductors can be classified in terms of a material parameter k > 0, known 
as the Ginzburg-Landau parameter. For the so called type-II superconductors (corre- 
sponding to k > \j\pi) there exist two critical magnetic field intensities, H Cl < H C2 , 
such that for \H S \ < H C1 the material is in the superconducting state and the magnetic 
field is excluded from the bulk of the sample except in thin boundary layers, while for 
\H S \ > H C2 the material behaves as a normal conductor, and the magnetic field penetrates 
it fully. For H Cl < \H S \ < H C2 a third state exists, known as "mixed state" as well as 
"vortex state" . The mixed state is characterized by a partial penetration of the magnetic 
field into the sample, which occurs, at a mesoscopic level, by means of thin filaments of 
normally conducting material carrying magnetic flux and circled by a vortex of supercon- 
ducting current. Type-I superconductors are those with H Cl = H C2 , so that the mixed 
state does not occur. 

Although from a physical point of view the most interesting description of the mixed 
state for type-II superconductors is given by the mesoscopic Ginzburg-Landau model, 
in designing magnets and other large-scale applications of superconducting materials, 
engineers use macroscopic models, involving averaged variables. One of the most reliable 



macroscopic models is the Bean's critical state model (see M). We refer to [13] for a 
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derivation of this model as a macroscopic version of the Ginzburg-Landau model under 
suitable assumptions. 

The basic idea of the Bean's phenomenological model for isotropic materials is that, 
because of physical limitations imposed by the material properties, the current density \J\ 
cannot exceed a critical value J c without destroying the superconducting phase. 

Moreover it is assumed that any electromotive force due to external field variations 
induces the maximum current density flow, according to the most effective way of shielding 
field variations. Hence the current density \J\ is forced to be J c in the part of the sample 
where the field is penetrated, while J = in the remaining part of the sample. 

At a mesoscopic level, the critical current density J c corresponds to the balance between 
a repulsive vortex- vortex interaction and attractive forces towards the pinning centers. At 
a macroscopic level, the electric field is zero when \ J\ < J c and abruptly rises to arbitrarily 
large values if J c is overrun. 

Since for isotropic materials it is well-established that the electric field and the current 
have the same direction (at least for slowly varying external fields), in the original Bean's 
model the Ohm's law is replaced by a vertical current-voltage law 

E\\J, \J\<J C , E = if \J\ < J c , 
which can be interpreted as the limit for p — > +oo in the power-law 



(1) E\\J, \E\=e c \±± 




(see, e.g., [10, ||, ||). Notice that the direction of the electric field is obtained by exploiting 
the isotropy of the material and cannot be obtained as a consequence of the power law 
approximation, which involves only the intensities of the fields. Hence this approach 
is not appropriate to deal with anisotropic materials (see anyhow || for some explicit 
computation of the electric field in the case of infinite slab geometry, and |8| in the case of 
a cylindrical body with elliptic section). On the other hand, many type-II superconducting 
materials are anisotropic. Moreover, even if the material is isotropic, the presence of flux- 
flow Hall effects has to be described in terms of an anisotropic resistivity (see, e.g., []Io|l). 

Concerning the constraints on the current, in the anisotropic case the Bean's law dictates 
that there exists a compact set A containing the origin as an interior point, such that J 
cannot lie outside A without destroying the superconducting state. Moreover J £ dA in 
the penetrated region, and J = elsewhere. 

Section ^ of this paper is devoted to the description of the Bean's law as a limit of a 
power-like law fulfilled by the dissipation. This approach allows us to determine a priori 
the direction of the electric field in terms of the direction of the current for anisotropic 
materials. 

In Section |] we deal with a variational statement of the critical state proposed in [Q], 
which takes the form of a quasistatic evolution of the penetrated magnetic field, obtained 
combining the finite-difference expression of Faraday's law with the Bean's law. We shall 
give a mathematical justification of the variational model as a limiting case of the power 
law model for dissipation. The result is proposed in terms of T-convergence of functionals, 
which is nowadays a classical tool in the mathematical methods for the material science 
(see, e.g., [||, [l6| and the references therein). 

In the last part of the paper we focus our attention to the special case of a long cylindrical 
anisotropic superconductor placed into a nonstationary, uniform axial magnetic field. The 
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parallel geometry enables us to make a two dimensional reduction of the problem, which 
can then be explicitly solved. 

The plan of this part of the paper is the following. In Section || we describe the two 
dimensional reduction of the problem in the case of parallel geometry. In Section ^ we 
introduce some technical tool needed for the analytical description of the fields. In Section 
|7] we find explicitly the solution to a general class of minimum problems with a gradient 
constraint and we determine the Euler equation solved by the optimal function coupled 
with its dual function. Moreover we find the explicit form of the dual function. In Section 
H the previous results are applied to the variational model for the critical state, and we 
find the quasistatic evolution of both the magnetic field and the dissipation inside the 
superconductor. Finally, a passage to the limit on the time layer gives the explicit form 
of the macroscopic electrodynamics. 

For what concerns the magnetic field, our result generalizes the one, valid for isotropic 
materials, obtained by Barrett and Prigozhin in Q with a different method based on a 
evolutionary variational inequality. Due to the fact that the magnetic field is explicitly 
known, we can compute the full penetration time in the case of monotonic external fields 
as well as we can depict the well known hysteresis phenomenon. On the other hand, the 
knowledge of the electric field allows us to give a detailed description of the evolution of 
the dissipation for cylinders of anisotropic materials with a general geometry of the cross 
section (see Figures [j] and |2|) . 



2. Notation and preliminaries 

For £ G R N , |£| will be the Euclidean norm, and (£, £') will denote the scalar product 
with £' G R . The symbol x will be used for the cross product of vectors in R 3 . Given 
a, b G R, a V b and a A b will denote respectively the maximum and the minimum of a and 
b. 

Given A C R N , we shall denote by Lip(A), C(A), C b (A) and C k (A), k G N the set of 
functions u: A — » R that are respectively Lipschitz continuous, continuous, bounded and 
continuous, and fc-times continuously differentiable in A. Moreover, C°°(A) will denote the 
set of functions of class C k (A) for every k G N, while C k,a (A) will be the set of functions 
of class C k (A) with Holder continuous fc-th partial derivatives with exponent a G [0,1]. 
Finally IP (A), W l ' p {A), Wq' p (A), and V(A,R d ), W ljp (A, R d ), Wq' p (A,R d ) , d > 1, will 
be the usual Lebesgue and Sobolev spaces of scalar or vectorial functions respectively. 

Let X be subset of R N , N > 2. We shall denote by dX its boundary, by intX its 
interior, and by X its closure. The characteristic function of X will be denoted by xx- 
The set X is said to be of class C k , k G N, if for every point xq G dX there exists a ball 
B = B r (xo) and a one-to-one mapping ifr: B —* D such that tp G C (B), ip -1 G C k (D), 
ip(B n X) C {x G R n ; x n > 0}, tp(B n dX) C {x G R"; x n = 0}. If the maps ip and 
are of class C°° or C k ' a (k G N, a G [0, 1]), then X is said to be of class C°° or C k ' a 
respectively. 

Let D C R^ be a compact convex set containing as an interior point and with 
boundary of class C 2 . The gauge function of the set D is the convex, positively 
1-homogeneous function defined by 



p D (0 = inf{t > 0; £ G tD} , £ G R 
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Since D is a compact set containing a neighborhood of 0, there exist two positive constants 
c\ < C2 such that 

(2) Ci |£|<Pd(£)<c 2 |£|, V?Gl w . 

The indicator function of the set D is defined by 

(O if£GL> 



(3) IniiJ 

Since D has a smoo 
plicitly computed, and 



Since D has a smooth boundary, the subgradient of the indicator function can be ex- 



{aDp D (£):a>0} if f G 3D 
(i) ^,(0-^0 if£££> 

{0} if £ G int D 



(see e.g. [22|, Section 23). 

In what follows, a family of objects, even if parameterized by a continuous parameter, 
will also be often called a "sequence" not to overburden notation. 

3. The physical setting 

Since we shall deal with a macroscopic model, the coarse-grained electrodynamics will 
be formulated in terms of 

(i) the flux density B within the sample, which is the average of the microscopic field 
intensity; 

(ii) the magnetic field H, which is assumed to be linearly connected to B by B = ^oH; 

(iii) the averaged current density J, linked to H by the Ampere's law J = curl H . 
Moreover, on neglecting finite size effects, we can assume that 

(iv) the magnetic source H s enters as a boundary condition for the flux density at the 
surface of the sample, requiring that n x H = n x H s on that surface, where n is 
the outward normal vector. 

Time variations of the external field H s induce in the sample an electric field E, ac- 
cording to Faraday's law. On the other hand, the electric field leads a current J which 
induces an internal magnetic field H, according to Ampere's law. Finally, we recall that 
div.ff = (Gauss' law). 

Summarizing, we are considering the so-called eddy current model for Maxwell equa- 
tions: 

-> dH 
curli? + ^o~q^ = (Faraday's law), 

J = curli? (Ampere's law), 

divi? = (Gauss' law), 

n x {H — H s ) = on the surface. 

In order to complete the physical setting, it remains to find an appropriate version of 
the constitutive law E(J). Our starting point is a reading of the macroscopic behaviour of 
the material in terms of the dissipation S = ^E(J), Jj (here 5 denotes the entropy of the 

system). Namely, the Bean's model dictates that if J is in the interior of an allowed region 
A, which is assumed to be a convex compact set of M 3 having the origin as an interior 
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point, then S vanishes as no electric field is generated in stationary condition, while when 
J touches the boundary <9A of the allowed region a huge dissipation occurs, destroying the 
superconducting phase. The fact that S = if J is an interior point of A, while S = +00 
if J A suggests that the Bean's law can be interpreted as the limit as p — > 00 of a power 
law for the dissipation (^E p (J), Jj, that is 

(5) (E P (J), J) = £ ( PA (J)) P . 

This approximation allows us to recover the direction of the electric field E(J) when 
J G <9A. Namely, we claim that, under the physically consistent assumption that 



(6) lim + EpitJ) = , VJG 



p3 



the relation (||) implies, for p > 1, 

(7) E p (J) = - ( PA ( J))^ 1 Dp A (J) . 

In order to prove @, we differentiate dH), obtaining 

(L»4(J), J) + E P (J) = c (pa( J)) P_1 DpA(f) ■ 
Hence, fixed J 7^ 0, the function v(t) = E p (tJ), t > 0, is a solution of the O.D.E. 

\v(l) = E p (J), 
where w = (pa(.J)^) P Dp/^(J). Then 



Ep (J) + —(t*>-l] 
p 



and, by (|), 



E P (J) - - (pa(J)Y D Pa (J) = lim [E p (J) + — (f - 1)] = 0. 

p v / t^o+ p 

As a consequence of (0), the direction of E p (J) is given by Dp A (J), and it does not 
depend on p. In conclusion, recalling (Q), the anisotropic version of the Bean's law is the 
following. 

(Bl) There exists a convex compact set A containing the origin as an interior point and 

such that J(x, t) G A for every x G t > 0. 
(B2) If we denote by I A the indicator function of the set A, the constitutive law E{ J) 

is given by E G dI A (J). 

4. A VARIATIONAL MODEL FOR THE MIXED STATE 

Let A C K 3 be the region occupied by the superconductor, and dA its surface. We 
assume that A is a bounded, open subset of M 3 , and that dA is of class C > 1 . Moreover, 
we assume that A has no enclosed cavity. In our model the presence of cavities is not 
allowed, since the boundary condition in such a cavity cannot be given in terms of the 
external magnetic field. (The magnetic field must be constant on the boundary of every 
enclosed cavity.) In order to avoid the presence of cavities, we assume that A has second 
Betti number (see e.g. [[0| for the definition of Betti numbers). 
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Let T > be fixed. For every n G N + , we take a partition P r 



{t?} 



k(n) 
i=Q 



of the 



interval [0, T], and we set 5tf = tf — tf_ l . The least action principle proposed in BJ is the 
following: starting from an initial field profile Hf(x) in A at the time layer i, and under 
a small change of the external drive, the new profile Hf +l {x) at the time layer i + 1 is the 
unique solution to the minimum problem 



(8) 



mm 



V - Hi\ 2 dx; V G Hs(t2 +1 ) + X 2 , cmlV G A a.e. 



where, for every p > 1, 



X p = [V G L P (A,R 3 ); curlT? G L P {A,R 3 ), divT? 



0, n x V = Oon 



In [GJ the trustworthiness of this variational principle was motivated in terms of nonequi- 
librium thermodynamical principles in analogy with the ohmic case. In this section we 
shall give a mathematical justification, in terms of T-convergence of functionals, of the 
variational model as a limiting case of the power law model for dissipation. 

Our starting point is the fact that, assuming that the power law (||) holds true, the 
discretized version of the Faraday's law 

curl^u, + -^r(i?f + i - &i) = 



(9) 



est: 



is the Euler equation for the functional 



F P {V) 



-PA (curl + A 
AP 



v-m 



dx . 



VeH s (t? +1 ) + x p 



where A := 



Mo 



2c8t 



- > 0. 



In other words, it can be easily checked that the magnetic field Hf +1 is the unique 
minimum point of F p in H s (tf, + X p . It is clear that, up to a translation of a constant 
vector, it is equivalent to minimize F p on X p by changing Hf with Hf — H s (tf +1 ). 

The following properties of the spaces X p will be useful in the sequel (see p5| , Thm. 2.2] 
and fH, Thm - 

Theorem 4.1. Let A C M 3 be a bounded open set of class C ' . Let 1 < p < oo. Then 
the following hold. 

(i) Every function V G X p belongs to W 1,P (A, M 3 ). Moreover, if p > 3, then V G 
C°' a (A) and there exists a constant C p > 0, depending only on p and A, such that 



V 



curl V 



- <C P 

C°' a (A) ~ p 

(ii) // A has second Betti number 0, then there exists a constant M p > 0, depending 
only on p and A, such that 



DV 



< M r 



curl V 



for every V G X p . 

Remark 4.2. From Theorem 4.1 we infer that X p is a (proper) subspace of W 1,P (A, 

Moreover, for p > 3 the map defined by X„ 3 V >— > curly is a norm in X p equivalent 

v 

to the standard W 1,p norm. 
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From now on we shall always assume p > 3. In this case, by Theorem |4.l| (i), a function 
V G X p is continuous in A, and the boundary condition rax V is understood to be pointwise 
fulfilled. On the other hand, the divergence-free requirement on V G X p is understood in 
the sense of distributions, i.e. 



f V-D$ = 0, V* eC^(A,] 

J A 



For our subsequent considerations, we need to define all the functionals F p on the same 
Banach space. For this reason, naming Uq = Hf — H s (tf +1 ), we set 



1 



2 



(10) = <( Ja P 

_+oo, otherwise in I? (A 

We also define the functional 

/ A (curlT/) + |y-[7o| 2 (ix, ifVeX 2 , 

otherwise in L 2 (A 

where I a is the indicator function of the convex set A, defined in pi). 



f(v) 




Lemma 4.3. The functionals F p , p > 3, and F are lower semicontinuous in the strong I? 
topology. Moreover, the functionals F p are equicoercive, i.e. there exists a constant a > 
such that 

F P (V) > ^\\V\\l -a VV G L 2 (A,R 3 ), Vp > 3 . 

Proof. Let 3 < p < oo, and let (V^) C L 2 (A, M 3 ) be a sequence converging to V in L 2 . We 
have to prove that 

F p (V)<limmfF p (V k ). 

k— >oo 

Without loss of generality we can assume that Vk G X p for every k and 
(12) lim F„(t4) = C < +oo . 

fe^oo 

We claim that V G X p . Namely, the divergence-free requirement is stable under strong 
L 2 convergence. Moreover, from (12) and @ we have that there exists a constant C\ > 



such that || curl Vfc ||p < C\ for every k G N. From Theorem |4.l| (i) we infer that the sequence 
(Vk) is equicontinuous and equibounded in A. Hence we can pass to a subsequence, that we 
do not relabel, converging uniformly to V in A, so that the boundary condition n x V = 
on dA holds. Finally, from Remark |4.2| , (Vk) converges to V weakly in W 1 ' P (A,R 3 ), and 
|| curll/ ||p < liminffc || curl Vfc||p, proving the claim. 

By the convexity of /?a and the convergence of curl Vk to curl V in the weak FP topology, 
we obtain that 

/ p A (curiy) p < liminf / p A (curlF fc ) p , 
J A fc^oo J A 

which implies the semicontinuity inequality F p (V) < C. 

The semicontinuity of F can be proved in a similar way, while the equicoerciveness of 
(F p ) easily follows choosing a = f ||C/q|||- D 
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Remark 4.4. For every p > 3, the functional F p is lower semicontinuous in L (A,~, 
coercive, and strictly convex in its effective domain X p . Hence it admits a unique minimizer 
Up, which belongs to X p . Similarly, the functional F admits a unique minimizer U £ Xi 
with pa (curl U) < 1. 

The mathematical justification of the Badfa and Lopez model is based on the fact 
that the sequence (U p ) of minimizers of (F p ) converges to the minimizer U of F. As is 
customary in material science, we are going to prove this statement using the technique of 
T-convergence (see |j) and the reference therein for examples of applications of the theory 
to different models). 

Definition 4.5. The sequence (F p ) T-converges to the functional F (with respect to the 
strong L 2 topology) if the following two conditions are satisfied. 

(i) For every V G L 2 (A, M 3 ) there exists a sequence (V p ) (called a recovering sequence) , 
converging to V in L 2 (A,M. 3 ), and such that 

limsupF p (Vp) < F(V) . 

p— ->oo 

(ii) For every sequence (V p ) C L 2 (A, M 3 ) converging strongly to V, we have 

liminf FJVp) > F(V). 

p— >oo 

Remark 4.6. Since the functionals F p are equicoercive (see Lemma |4.3| ), the convergence 
of the minimizers (U p ) to U will follow if we prove that (F p ) T-converges to F and that 
(Up) is a precompact sequence in L 2 (A, IR 3 ) (see e.g. |J). 

The main tool needed in the proof of the r-convergence of (F p ) to F is the following 
lemma. 



Lemma 4.7. Let (V p ) G L (A,M. ) be a sequence converging to V in L (A, R ). // 

liminf F V (V V ) < +00 

p— >oo ^ 

then pa (curl V) < 1 a. e. in A. 

Proof. Without loss of generality we can assume that 

liminf FJV V ) = lim F P (V P ) = L< +00, 

p— >oo ^ r p— >oo r r 

curl V p e L P (A, M 3 ) , F p (V p ) < L x Vp > 3. 

We are going to prove that the sequence (V p ) is bounded in W l,q (A, M 3 ) for every q > 3. 
In view of Theorem 4.1, it is enough show that for every given q > 3 there exists a constant 
C > such that || curlV^,|| < C for every p > q. 

Given 3 < q < p, by (||) and Holder's inequality we have 

/ I curl V p \ q dx < -g [ pa (curl V p ) q dx 
J A q J A 

4(^) s/ M/> a(ci,ri ^r- 



(13) 

On the other hand we get 



- p A (cmlV p ) p dx < Fp(V p ) < L x . 
A P 
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Collecting all the previous estimates, for every p > q we obtain 
(14) /. | curled* < ^ f^V /P < ^ exp 



/a 1 pi " 4 V \A\ J - 4 ~*\e\A\ 

Then the sequence (V p ) is bounded in W 1,q (A, K 3 ), and hence (V p ) converges to V in the 
weak topology of W 1,9 (A,R 3 ). 

This fact and the convexity of the function p& imply that 

(15) / pA(curl V) dx < liminf / p/^(cxii\V p ) dx 

Jb p^°° Jb 

for every open set B C A (see e.g. (O, Thm. 4.2.1]). On the other hand, by Holder's 
inequality 

J PA (cuAV p ) dx < (^J PA (cmlV p ) p dx^j * \B\^/ p 



< V llv Q J PAicuvlVpfdx^ P {B^-V/p . 
Then, using (0) and (H) we get / /9a (curl V) < \B\. Since /^(curlVj G L 1 {A) 1 we 



conclude that 



PA (curl V(x)) = lim / pa (curl V) dx < 1 

r->0+ \BJx)\ JbJx) 



+0+ |S r (x)| JB r {x) 

for a.e. x G A. □ 
Theorem 4.8. T/ie sequence (F p ) I '-converges to F . 

Proof, (i) Let V G L 2 (^4,M 3 ). It is not restrictive to assume that F(V) < +oo, so that 
V G X2 and pA(curiy) < 1 a.e. Then it is enough to choose V p = V as recovery sequence, 
(ii) Let V p G L 2 (A,R 3 ) be functions converging to V in L 2 (A, R 3 ), such that 

liminf -K, (14) < +00. 

By Lemma |4.7| we have that pa (curly) < 1 a.e. in A, and hence 

F(V) = I \V-U Q \ 2 dx = lim / \V P - U \ 2 dx < liminf F P (K) , 
concluding the proof. □ 

Notice that the T-convergence result is obtained under the sole assumption that Uq G 
L 2 (A,M. 3 ). Actually in our model U = Hf - H s (tf +1 ) so that U e X 2 + c , c G M 3 , and 
Pa(cu.t\Uo) < 1. This additional regularity is needed in order to recover the convergence 
of the minimizers. 

Theorem 4.9. Let Uq G c*o + X2, c*o G M 3 , with pa(cut1Uq) < 1 a.e. in A. Let U p be 

the unique minimizer of F p , p > 3, and let U be the unique minimizer of F. Then (U p ) 
converges to U in the weak topology of W 1,q (A, M. 3 ) for every q > 3. 

Proof. Reasoning as in (jl^), for 3 < q < p we have 

(16) j A \cvLT\U p \Ux<^ (jjyj q/P \A\ y A ^ P A(cmlU p ) p d x y P 
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On the other hand, by the minimality of U p we get 

/ - pa(cut\U p ) p dx < F p (U p ) < F P (U - c ) < (1 + A |c | 2 )|^| 
Jap 



so that 



/ \ cm\U p \Ux< ^[p(l + A|c | 2 )] g / p |^| <\\A\ exp[q(l + A |c | 2 )/e] . 
J A C\ C\ 

Then the sequence (Up) is bounded in W 1,q (A, R 3 ). In particular, it is precompact in 
L 2 (,4,R 3 ). Then, by Remark [D we obtain the convergence of (U p ) to U in L 2 (A,\ 



Finally, notice that the previous estimate implies that (U p ) converges to U in the weak 
topology of W 1,q for every q > 3. □ 

5. Cylindrical superconductors 

In what follows we shall consider a long cylindrical type-II superconductor occupying 
the region fixR, with a simply connected cross section !!cK 2 , placed into a nonstationary 
uniform axial magnetic field H s (t) = (0,0, H s (t)). 

Notice that, due to the parallel geometry, the problem admits a dimensional reduction. 
Namely, we can assume that there exists a function u depending only on x = (x\,X2) S 
VL such that H = (0,0, H s (t) + u(x,t)), so that curli? = (u X2 (x, t), — u Xl (x, t), 0) and 
div H = 0. In particular the current density J is parallel to the cross section plane and 
its projection on this plane is (Ji, J2) = (Du) ± , where (Du) ± is the rotation of — | of 
Du. Hence, under the additional hypothesis that the allowed region A is symmetric with 
respect to the z = plane, the constitutive law becomes Du € K, where K is the rotation 
of the section z = of the set A, and E = (Ei,E 2 , 0) with (Ei,E 2 ) £ 81k (Du). 

In this case, if c" = H s (tf), and hf(x) is the magnetic field intensity at time layer i, 
given by H s (tf) + u(tf,x), the variational formulation is 

(17) min jy (w - hf) 2 dx; w G c™ +1 + Wq'^O), Dw 6 K a.e. in n\ . 

In the following sections we shall give an explicit representation of the unique minimizer 
hf +1 and we shall show that the previous quasistatic evolution converges to a function 
h(x,t) (also explicit), which determines the evolution of the internal magnetic field. 

Further to the dimensional reduction and a rescaling, we obtain that hf +l — c" +1 is the 
unique minimum point of the functional 

J(v) = f \l K (Dv) + (v- u) 2 ] dx, v G Wo A (n), 
where u = hf — c" +1 . 

Although in this case we deal with an unbounded superconductor A = £1 x R, the above 



reduction allows us to prove, using the same arguments as in Section [| (see also [18], 
Section 2), a T-convergence result, and, consequently, the convergence of the minimizers. 
Even if the case iV = 2 is the sole meaningful physical situation, the following results hold 
for every dimension N > 2. 

Theorem 5.1. Let Q, C M. N be a bounded open set, let K C ~M> N be a nonempty, compact, 
convex set with G hit if, and let px be the gauge function of K. Let G p , p > 1, and G 
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be the junctionals defined respectively by 
G,Ar) 




p K (Dv) p dx, veW^ p (fl) 

otherwise in L 1 (J7) 



and 



jO, v G W 1,oc (fi), p{Dv) G if 
1 +00 otherwise in L 1 (J1) . 

Then (G p ) T -converges to G in L r (Q) for every r > 1. As a consequence, for a given 
u G L 2 (J7) i/ie functionals 



G(v) 



J p (v) 



n 

+00 



1 , _ 2 

- pk(Dv) p dx + (v — u) 
P 



dx, v G W 1,p (ft) 

otherwise in L 2 (£l) 



r -converge to the functional 



J(v) 




I K {Dv) + {v-u) 2 dx, v G Wq'^O) 



i,i/ 



otherwise in L 2 (£2) , 

in L 2 (Q). If in addition u G Co + VFq' 00 ^), Co G R, imi/i p{Du) < 1 a. e. in $7, i/ien £/ie 
minimizers u p of J p converge to the unique minimizer u of the functional J in the weak 
topology 0/ Wg' 9 (fi) f or every q > 1. 

6. The Minkowski distance function 
Here and hereafter we shall deal with 
(18) U C R N nonempty, bounded, open connected set of class C 2 , 

and 

K C M N nonempty, compact, convex set, G hit if, 
dK of class C 2 with strictly positive principal curvatures . 
The polar set of K is defined by 

K° = {p G R N ; (p, x) < 1 Vx G K} . 
We recall that, if K satisfies ©, then K° also satisfies ©, and if 00 = (if ) 



(19) 



K (see 
fH, Thm. 1.6.1]). 

Since K will be kept fixed, from now on we shall use for the gauge functions the notation 
p = PK and p° = p K o 

It can be readily seen that the gauge function p° of the polar set K° coincides with the 
support function of the set K. As a consequence, we have that 

(20) p(Dp°(0) = l V£<=R N \{0}. 
Let Lip p (f2) be the set of functions defined by 

Lip p (f2) := {u G Lip(H) : Du G K a.e. in Q}. 
We recall that u G Lip p (Q) if and only if 

(21) u(x) - u(y) < p°{x - y) , 
for every x, y G 0, joined by a segment contained in Q. 



12 



G. CRASTA AND A. M ALUS A 



The main tools needed in the following sections are the distances from the boundary of 
Q associated to the Minkowski structures induced by the gauge function of K° and — K° 
respectively. 

Definition 6.1. The Minkowski distance from the boundary of Q is 

(22) d n {x) = inf p°(x - y), xgH. 

Similarly, we define d^(x) as the distance from the boundary induced by the gauge of 
-K°. 

Notice that the function coincides with dn only if K° is symmetric with respect to 
the origin. In the remaining part of this section we illustrate some features of the function 
dn- The analogous for d^ can be obtained upon observing that — K° is the polar set of 
—K. 

Since dQ is a compact subset of W N and p° is a continuous function, the infimum in the 
definition of dn is achieved. We shall denote by II (x) the set of projections of x in dQ, 
that is 

(23) II(x) = {y G dn- d n (x) = p°(x - y)}, x G O. 

Definition 6.2. We say that x G is a regular point of 0, if II(x) is a singleton. We say 
that x G f2 is a singular point of Q if x is not a regular point. We denote by S C the 
set of all singular points of f2. 

Remark 6.3. It is well known that dn is a Lipschitz function in S7, and that do, is dif- 
ferentiable at x G fi if and only if x is a regular point. In addition, if II(x) = {y} 
then Ddn(x) = Dp°(x — y), and hence, by 
|, 0, (H); see also Theorem |J(i) below). 

Remark 6.4. We recall that do, (resp. d^) is the unique viscosity solution of the Hamilton- 
Jacobi equation p(Du) = 1 (resp. —p(Du) = —1) in £1, with boundary condition u = on 
dd. As a consequence of the maximality property of the viscosity solutions, we have that 

— (Iq < u < dn , Vn G Lip p (Sl), n = on 9fi 

(see [^l}). Moreover it can be easily checked that 

u > — , Vn G Lip p (Sl), u > on , 

u < dn , Vn G Lip p (0), n < on , 
(see e.g. 0], Theorem 5.9). 

In the following theorem we collect all the results proved in |l4j that are relevant for 
the subsequent analysis. 

Theorem 6.5. Assume that Cl and K satisfy respectively ( |i~8| ) and (|l^). Then the following 
hold. 

(i) E C f2, and the Lebesgue measure of £ is zero. 

(ii) iei x G SI and y G II(x). Then dn is differentiable at every z along the seg- 
ment jointing y to x (without endpoints), and Ddn(z) = ^jj—yy, where v{y) is the 
(Euclidean) inward normal unit vector to dVL at y. 

(iii) The function dn is of class C 2 in Q, \ £. 

Proof. See Remark 4.16, Corollary 6.9, Lemma 4.3 and Theorem 6.10 in [^]. □ 



2C), Ddn G K almost everywhere in Q (see 
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At any point y £ d£l there is a unique inward "normal" direction p(y) with the properties 
II(y + tp(y)) = {y} and dn(y + tp{y)) = t for t > small enough (see |[4[ Remark 4.5]). 
More precisely, these properties hold true for p(y) = Dp(v{y)) and for every t £ [0,l(y)), 
where l(y) is defined by 

(24) l(y) = mm{t > 0; y + tDp(Dd n {x)) £ £} 

(see Jl4], Propositions 4.4 and 4.8] and [2C, Lemma 2.2]). From Theorem |6.5| (ii) and the 
positive 0-homogeneity of Dp it is plain that Dp(v{y)) = Dp{Dd,Q j {y)). 

From Theorem [0](iii), the function do, is of class C 2 on <9$7. Then we can define the 
function 

(25) W(y) = -D 2 p(Ddn(y))D 2 d n (y), y£dn. 

For any y £ <9$7 let T y denote the tangent space to d£l at y. It can be proved that for 
every v £ T y , one has W(y)v £ T y . Hence, we can define the map 

(26) W(y):T y ^T y , W(y)w = W(y)w, 

that can be identified with a linear application from IR^™ 1 to l^^ 1 . 

Although the matrix W(y) is not in general symmetric, its eigenvalues are real numbers, 



and so its eigenvectors are real (see [14], Remark 5.3). The eigenvalues of W(y) have an 
important geometric interpretation. 

Definition 6.6. Let y £ d£l. The anisotropic curvatures of dfl at y, with respect to the 
Minkowski norm dn, are the eigenvalues K\(y) < ■ ■ ■ < £ n _i(y) of W(y). 

In what follows we shall denote by II - , k~ , j = 1, . . . , n — 1, X - and l~ , respectively 
the projection, the anisotropic curvatures and the singular set associated to d^, and the 
normal distance to the cut locus £ _ . 



7. A MINIMUM PROBLEM 

Fixed u £ Lip p (f2), we want to minimize the functional 

(27) J(v) = f [I K (Dv) + (v-u) 2 ]dx, v£W^ 1 (n). 

Let us consider the partition $7 = £l + U U with 

S7° = {-d^ <u< dn} , tt + = {u > cfo} , = {u < -d^} , 
and define the function 

!d n (x) x£fl+ 
-d^{x) x £ U~ 
u(x) x £ U° . 

The aim of this section is to show that u is the unique minimizer of J in VFq' 1 ^). 

Remark 7.1. If x £ , then u{x) = tfo(x) = — d^(x) and hence x £ dQ andu(x) = 0. 

As a consequence, we have u(x) = min((max(u(x), — d^(x)), cfo(x))). This implies that 
u £ Lip p (f2), u = on dQ, p(Du) = 1 a.e. in + U and u is a viscosity solution of 
p{Du) = 1 in Q + and of —p(Du) = — 1 in Q~ . 
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(30) 



Let us define the following two subsets of d£l: 

(29) r+ = {y G 8Q: u(y) > 0}, T~ = {y G 3ft: u{y) < 0} , 

and the functions 

A(y) = sup{s G [0, Z(y)) : u(y + tDp{y{y))) > t V t G [0, *)} , 

A" (y) = sup{s G [0, r (y)) : «(y + tDp. {u{y))) < -t V t G [0, s)} . 

Recall that Z(y) and Z _ (y) are the normal distance to the cut locus, corresponding to da 
and d^ respectively, that is the length of the "normal" ray starting from y and ending on 

S(see(H)). 

Remark 7.2. It can be easily checked that A and A~ are lower semicontinuous function in 
<9fi. In general we cannot expect the continuity of these functions, as it is shown by the 
following example. Let ft = Bi(0) C R 2 , K = £?i(0) C R 2 , and define u(r cos 0, r sin 6) = 
1 — rcos6> for < r < 1, G [0, 27r). We have |Du(r cos 6*,r sin #)| 2 = cos 2 8 + sin 2 6 = 1 



(0 < r < 1), so that u G Lip p (ft). Moreover, if y 
lower semicontinuous function given by 



(cose, sine) G <9ft, 9 G [0,2tt), A is the 



0, 6> = 0, 

1, #G(0,2vr). 



The following simple lemma provides a characterization of ft + and ft in terms of the 
rays starting from points of r + and T~ respectively. 

Lemma 7.3. Let y G T + be such that A(y) < l(y). Then y + tDp{v(y)) G ft for every t G 
[A(y),Z(y)]. Analogously, lety G T~ be such that \~ (y) < l~(y). Theny+tDp_{y(y)) G ft 
/or every t G [A~(y), Z~(y)]. ^4s a consequence, we have 

ft+ \ £+ = {y + ^p(i/(y)) : y G T+, t G (0, A(y))}, 

A" \ E" = {y + tDp-(y(y)) : y € T" i G (0, A"(y))}. 



(31) 



Proof. Let y G T + and assume that A(y) < Z(y). From the very definition of A(y) and the 
continuity of u, we have that u(y + A(y) Dp{v(y))) = A(y). For every t G [A(y), l(y)] we 
have that 

dn(y + tDp{v{y))) = t = u{y + A(y) Dp(u(y))) + t - A(y) > u(y + i Dp(u{y))) , 

where the last inequality follows from the assumption u G Lip p (ft) (see (|2l|)), and from 
the fact that p°((t - \)Dp{v{y))) = t- A (see @)). The inequality u(y + tDp(u(y))) > 
—d^(y + t Dp(u(y))) follows easily from the fact that u(y) > and u G Lip p (ft) (see 

□ 



Remark |6.4| ). The case y G T can be handled in a similar way. 

The following change of variable formula can be proved as in [14, Theorem 7.1]. 
Theorem 7.4. Let <£, <I>~ : <9ft x IR ^ R N be the maps defined respectively by 

<S>{y,t) = y + tDp{u{y)), $~ (y,t) = y + t Dp_{v{y)), {y,t)edQxR. 



Then 



/ h(x) dx = 


f pWv)) 


11 


/«+ 






/ h(x) dx = 


j r _P-{v{y)) 


In- 





n-1 



A" (y) 



i=l 



n-1 



dH n -\y) 

m—l i 



dU n -\y) 
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for every h G L 1 (r2). 

In what follows we will set 



n— 1 



MS) = n i 



1 — tki(x) 



n-l 



, =1 . da(x)ki(x) ' Mx (t) 2 ! _ d n {x)RT( x ) 
From Lemma 7.3 in |l4| we deduce that there exists a positive constant Mq such that 
(32) < M x (t) <M ViGfi\ E, d u (x) < t < l(y), 

where LI(x) = {y}. An analogous estimate holds for M~(t). 



1 - tie* (x) 



Using Theorem 7A and (|31|), we are able to show that the function u defined in (|28[) , 
coupled with an explicit function v, solves a system of PDEs of Monge-Kantorowich type. 
A posteriori this system will be understood as a first order necessary condition fulfilled 
by the minimizer of J. 

Lemma 7.5. Let u G Lip (Q) be the function defined in (28), and let v. £1 —* M be the 
function defined by 
(33) 



v(x) 



m 

\ 

dn{x) 
A- (y) 



u(y + tDp(Ddn(x)))-t]M x (t)dt, x G n+ \ E 

[u{y + tDp-(Dd n (x))) +t]M~(t)dt, x G Q~ \ £~ 



n(x) = {y} 
U~(x) = {y} 



otherwise . 

Then v G Cf,(f2) ; w > 0, and £/ze pair (u, v) solves the system of PDEs 

— div(v Dp(Du)) =u — u in ft, ( distributional) 

p(Du) < 1 a.e. in $7, 

p(Du) = 1 a.e. in {v > 0}, 

u = in 917 . 



(34) 



Proof. Let us extend the functions A and A - to Q, + \ E + and fi _ \ E~ respectively by 
setting A(x) := A(y) when LI(x) = {y} and A~(x) := A~(y) when n _ (x) = {y}. Since 
A(x) > dn(x) for every x G Sl + \ E + and A~(x) > d^(x) for every x G \ £~, it is plain 
that v > 0. 

Let us prove that v is a continuous function. It is convenient to rewrite v in f2 + \ £ + 
in the following way: 

r«v) , _+ 

(35) v(x)= [t\/u(y + tDp(Ddn(x)))-t]M x (t)dt, x G fi \ £ , LI(x) = {y}. 

0) 

Since the maps cZq, DdQ, u, I, and x i— > LI(x) are continuous in £l + \ £ + , it follows that 
also v is continuous in f2 + \ S + . 

If xo G ri + H £ , then A(xo) = ^(xo) = dn(xo), and the map A is continuous at xo- 
Namely, since xo G £ , we have that l(xo) = dn(xo). Moreover, since xo G f2 + , we have 
that H(xo) > dn(xo), so that u > cfo on every segment [y, xo] C with y G 50. Thus 
A(xo) > dn(xo). Finally, from the lower semicontinuity of A and the fact that A(x) < l(x) 
for every x G 0, we conclude that 

dn(xo) < A(xq) < liminf A(x) < limsupA(x) < lim l(x) = 1(xq) = dn(xo), 
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so that \(x ) = du(xo) = lim x -> xo X(x). 



From the very definition of v and ( |32|) we have that 

< v(x) < C[X(x) - d n (x)] Viefi + \ S + . 
We remark that the same estimate trivially holds also on fi + n S + since v vanishes. If 



xq G VL + n £ + , we have that 

limsupt^x) < limsupC[A(x) — dn(x)] = 0, 

X— -rXQ X— >XQ 

hence lim a ._ +xo v(x) = = v(xq). 

Up to now we have proved that v is continuous in fi + . A similar argument shows that 
v is continuous in Q~ . Since v vanishes on 0° and dfl + fl dQ~ nH = (see Remark 7.1), 
it remains to prove that v is continuous onflfl d£l + and Q, fl . 

Let xq G H <9f2 + . Since xo G 0°, we have that u(xo) = 0. Moreover, by definition of 
Q + , u(xo) = dfi(xo). Since u G Lip p (f2), by ( f2~I| ) and (||), given e > we have that 

u{x) < dn(x) + 2c2e Vx G fin I? e (xo). 

Moreover, if x G (fi + D B t (xo)) \ S + , we have that u < dn + 2c2e along the segment joining 
x to its cut point. As a consequence, the integrand in (^) is bounded from above by 
2c2£ Mo, so that 

0<v{x)<Ce \/x G fl B e (xo)) \ S + . 
Since -u = in U S + , this inequality clearly implies that 





lim v(x) 

X—>XQ 



V(Xq) 



i.e. v is continuous at xo. The continuity in £1 n <9fi can be proved in a similar way. 
Let us prove that the pair (u,v) is a solution to (Q). Upon observing that {v > 0} 
U fi~, by Remark 7.1 we have only to prove that (u,v) solves 

(II — u)ipdx , 



v(Dp(Du), Dip) dx 
for all ip G Cg°(fi). We have that 



(u — u)ip dx 



(u — u)ip dx + (u — u)ip dx 



Q- 



By using the change of variables stated in Theorem [7.4| , we get 

r \- n-l 



n- 



(u — u)f 



dU 



n-l 



An integration by parts leads to 

rX 



« - i 



A" 



n-l 



1=1 



Since u($ (y, s)) : 



s for every s£ [0, A (j/)], we have that 

X-(y) n-l 



(u^-(y,s))-u^-(y,s))) \[ \ <!■ 

i=i 1 



*«i (y) 
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for every y £ T and t € [0, A (y)). Then we obtain 

/ (u-u)ipdx = - [ p-(v) [ X (D<p($-), Dp-{v))v{$) T\(l-tKi)dt 
Jn- Jr- [Jo " 

Finally, again from Theorem |7.4| and the fact that 

Dp.(u(y)) = -Dp(-u(y)) = -Dp(-Dd n (y)) = -Dp(Du(y)) , 

we get 

(u — u)(pdx = I v (Dp(Du), Dip) dx , 



dU 



n-l 



n- 



for every ip £ C£°(Q,). In the same way, using the change of variables in we get 

(u — u)ip = v (Dp(Du), Dip) , 



for every (p £ C£°(Q), concluding the proof. □ 
Now the absolute minimality of u follows easily. 



Theorem 7.6. The function u defined in (28) is the unique minimizer of J in W 1 ' 1 (Q). 



Proof. Let v be the function defined in fl33|). Since v > 0, for almost every x G + U 0~ 
the vector 2v(x)Dp(Du(x)) belongs to the subdifferential of Ik at Du(x) (see (H)). On 
the other hand, v = in 0°. Then we have 

I K {Dw(x)) - I K (Du{x)) > 2v(x) (Dp(Du(x)), Dw{x) - Du{x)) a.e. in O, 
for every w G VFq' (O). Hence, using Lemma |7.5| we obtain that for every w S W ' (O) 

J(w)-J(u)> [ 2v (Dp(Du), Dw - Du) + / (w-u) 2 + [ (u-u) 2 
Jn Jn Jn 

(36) = ( \2{w-u){u-u) + {w-u) 2 -{u-u) 2 } 

l 

! 

which implies that it is the unique minimizer of J in Wq^^VI). □ 

8. The quasistatic evolution 

Let T > be fixed. For every n £ N + , we take a partition P ra = {£"}*J.q of the interval 
[0, T], and we set 5if = tf - tf_ v 

The starting configuration of the superconductor at time t = is given by h(x, 0) = 
/io( 2; ) ^ Lip- (n), with /io = H s (0) on Using the variational formulation §L7\ ) and 
Theorem |7.6| , we have that the internal magnetic field hf +1 (x) = h(x,t™ +1 ) is given by 

(37) K? +1 (x) = mm[max[h?(x),<? +1 - d^(x)],cf +1 + d n (x)\ . 

Since the model is valid only when the external field varies slowly in time, it is not 
restrictive to assume that H s is a continuous function obtained by the concatenation of 
a finite number of monotone functions of class C 1 . Hence from now on we shall consider 
only the case 

(38) H s is a monotone function of class C 1 ([0, T]). 
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In this case formula ( |37|) simplifies. Namely, if we assume that H s is a nondecreasing 
function we have that cf < cf +1 . In particular, at the first time step we have 

hi(x) = min[max[/io(x), c™ — cIq(x)], c\ + dn(x)] , 

and Hq(x) — c" G LipJO) with h Q (x) — c™ = Cq — c" < on 90. Hence /io(x) — c™ < G?n(a;) 
in f2 (see Remark 6.4), which implies 

hi(x) = max[/io(x), — <i^(x)] . 
At the second time step we get 

h%(x) = min[max[/i"(x), C2 - d^(x)],c^ + dn(x)] . 
Since hi(x) — c 2 = max[h Q (x) — c%, — d^(x) + c" — c^], with c" — cJj < 0, we have that 

^2 ( x ) = max[/io(x), C2 — dft(x)] . 

In general we get 

(39) h-2(x) = max[/io(x), c" — c^(x)] . (nondecreasing external field) 

With a similar argument, we obtain that if H s is nonincreasing in time, the i— th time step 
leads to 

(40) ^r( x ) = min[/io(x), + dn(x)] . (nonincreasing external field) 

Now we easily obtain the evolution of the internal magnetic field, taking the limit as n 
goes to +00 of the quasistatic approximation 

h n (x,t) = h?(x), te [t?,t? +1 ). 

Theorem 8.1. If H s is a nondecreasing function in [0, T], then the sequence (h n (x,t)) 
converges uniformly in £1 x [0, T] to 

h(x,t) = max[/io(x), — d^ (x) + H s (t)] . 

If H s is a nonincreasing function in [0, T], i/ien i/ie sequence (h n (x,t)) converges uniformly 
infix [0, T] to 

/i(x,t) = min[/i (x), <fo(x) + #s(*)] • 

Proof. Assume that H s is nondecreasing (the other case being similar). Given t G [0, T], 
and u 6 N, we have that t G [if ,tf +1 ), and 

< h(x,t) - h n (x,t) < H s (t) - H s {tf) < H s {tf +l ) - H s (t?) < 6t n maxH' s (t) , 

where 5t n is the size of the partition P n . Hence the uniform convergence in Vt x [0,T] 
follows. □ 



Remark 8.2. Theorem 8.1 shows that the evolution obtained by the direct optimization 
method at discrete times proposed by Badi'a and Lopez in (l| coincides with the one 
obtained in |(| by Barrett and Prigozhin as a solution of an evolutionary variational 
inequality. 

For what concerns the electric field induced in the superconductor, we already know that 
E = (Ei,E 2 , 0) with (Ei,E 2 ) G dI K (Dh). This means that E = whenever Dh G intK, 
and there exists w(x,t) > such that (Ei,E 2 ) = w(x,t)Dp(Dh(x,t)) for almost every 
x G ft and for every t G [0, T}. Our approach allows us to compute explicitly the quasistatic 
evolution of the dissipated power w(x,t) as the dual function appearing in the necessary 
conditions stated in Lemma fill. 
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Namely, from Lemma 7.5 , for every n G N and for every i = 1, . . . , k(n), there exists 
vf £ Cb(£l) such that 



(41) -div I^Dp(DK, 



(42) 



Jtf rK tJ J St? ' 

in the sense of distributions in Q. For our convenience, we also set Vq = 0. The function 
vf is unique, as stated in the following result. 

Proposition 8.3. Let v £ Cb(Q) be a function such that 

{ — div (v Dp(Dhf)) = hf_ 1 — hf , in O (distributional), 
p(Dhf) = 1 a.e. in {v > 0}. 

If cf > cf_ 1} then v = vf , where 

1/ (h?-ht 1 )(y + sDp^(u(y)))M-(s)ds, xen\E ,IL-(x) = {y}, 

V { (x) = < Jd-(x) 

[o ieE" 

whereas, if cf < cf_ x , then v = vf , where 

( rKv) _ 

1/ (/>?. 1 -ft?)(j/ + 8iJ/»Ky)RW& ) I en\E,n(i) = { y } 1 

[o ieE. 

Proof. Assume that cf_ x — cf > (the other case being similar), and let v S Cb(fi) be a 
solution of (|42"1). By construction we have that the function / := hf_ x — hf is bounded, 
continuous and nonnegative in Q, and hf = da in supp(/). Using the same arguments of 
Lemma 6.3 in [15], we can prove that the pair {d^,v) satisfies 

-dw(vDp(Dd n )) = f, 

in the sense of distributions in fi. Now we can apply Propositions 6.5 and 6.7 in [15] in 
order to conclude that v = vf . □ 



Formula (Ell) can be understood as a discretized version of Faraday's law, so that the 



V 

functions wf = -A are the steps of the quasistatic evolution of the dissipated power. If 



we set 

(43) w n (x,t)=wf(x), g n (x,t)= ^ >— lii , t£[tf,tf +1 ) 



and w n (t) = w n (-,t), h n (t) = h n (-,t), g n (t) = g n (-,t), formula ( ]4T|) can be rewritten as 
(44) - div (w n (t) Dp{Dh n {t))) = g n {t). 

Theorem 8.4. If H s is a nondecreasing function in [0, T], then the sequence (w n (x,t)) 
converges pointwise in x [0, T] to 

dH s r l ~(y) 

W{x,t) = -K-{t) X{ho(y+sDp_(u(y)))<H 3 (t)-s}M x {s)ds. 

IfH s is a nonincreasing function in [0, T], then the sequence (w n (x,t)) converges pointwise 
infix [0, T] to 

dH s f l (y) 

W ( X > *) = STW / X{h (y+sDp(u(y)))>H s (t)+s}M x {s) ds . 

Ol Jda{x) 



20 G. CRASTA AND A. MALUSA 

Moreover, in both cases, the sequence (w n ) converges to w in the strong topology of L P (Q,), 
p>l, uniformly in [0, T\. 

Proof. We compute the limit of the sequence (w n ) only in the case of nondecreasing 
external field, the other case being similar. 

We can assume, without loss of generality, that for every n £ N the partition P n+l is a 
refinement of the partition P n . Hence, for every t £ [0, T] and for every n £ N there exists 
i = i(n) £ N such that t £ [i" (n) , £"( n \ +1 )- Moreover, since H s £ C 1 ([0,T]), we obtain 

(45) 7r^-{t) = lirn — ^ - — s ^ i uniformly in [O.Tl. 

(jt n-*oo dtl 



By Proposition |8.3| we get 



rl~(v) 

w n (x,t) = / g n (y + sDp-(u(y)),t)M-(s)ds, 
Jd n (x) 



where g n is the function defined in (|43|). Moreover we have 

K_ x -h? = (il s (d) - ^( t i l ))X{ ho <H s (t"_ 1 )-d n } 

+ (h - H s (t?) + rfn)X{ //s ( t n_ l) _ d - < ^ 0<//s(t n ) _ d -} , 



(46) 
and 



> (h - H a ($) + rfQ)X{ ffs(t? _ i )_ d - </l0<ifs(i?) _ d -} 
> - H s (t?)) X{ Hs(t? i) - dn<ho< H s (t?)~d n } • 

Collecting the previous information, we get 

f) J-f 

and lmin^+oo ||<7 n (i) — g(t)\\i P ^ = uniformly in [0,T], for every p > 1. In addition, 
from (HH) and fl46|), for every regular point x of and setting {y} = H(x), we get 



(48) lim g n (y + sDp_(u(y)),t) = g(y + sD P -(u(y)),t) a.e. s £ [d n (x), I (y)) . 

n—*+co 

Finally, by (^), (^) and the Dominated Convergence Theorem, we conclude that 

lim w n (x,t) = 
n— >+oo 

(49) dH s . . r l (y) . , , . . 

S-W /_ X{^ (3/+sDp_(i/(y)))<i? s (t)-s}-Mi; (*) : = ^(x, i) . 

The last part of the theorem follows from the fact that 

(50) \\w n (t) - w(t)\\ LP(n) < C \\g n (t) - g(t)\\ LP(n) 

for some positive constant Co- Namely, for a given t £ [0,T] let us define 
<p n (x, s) := g n (y + sDp_(v(y)), t) - g(y + sDp^(v(y)),t), 

for x £ f2 \ X, n(x) = {y}, and s G [0, l~(y)]. We remark that, if x £ f2 \ £ and y is the 
projection of x, then y> n (x, •) = <^ n (y, •). 

By the very definition of w n , w and M x , Holder's inequality and ( |3^ ) we have that 

r ( r^~( x ) n—i -. — / \ \ 

i-w - wwr„ (11) < o/ n / i^(„.)r n ..^^ * - ' 
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where I (x) := I (y) when IT (x) = {y} and C := Mg 1 (maxd^) 1 ^'. Using the change 
of variables theorem (see jl4|, Thm. 7.1]) we have that 







I = C[ p_(v(y)) 




J an 





n-l 



\^ n (y,s)\ p n i 



n-l 



i=l 



c p-Hy)) 

Jan 



I \ V n (y,sWl[(l-sk-(y))dsj da 

<C'J p-{u{y)) J \<p n (y, S )|f J] (1 " *«T(y)) ^"^(y), 



n-l 



1 - a«j (y) 
<™~(y) 



ds 



dn n 'Hy) 



where C := C maxd^. Using again the change of variables theorem we finally get 

/<C / |s n (x,t) -</(*, t)|*dx, 



hence 



follows. 



□ 



Remark 8.5. Notice that if d^(x) > A" (y), then (x) = 0, while /if (x) = c" — d^(x) if 
d n ( x ) ^ Hence we have that w n (x,t)Dp(Dh n (x,t)) = w n (x,t)Dp(-Dd^(x)) a.e. 

in f2 and for every t G [0, T]. A passage to the limit in ([lij) leads to 

(51) -divMa?, t)£>p(Z>h(s,t))) = — ^ X^^)^} ' 

in the sense of distributions in £1 for every t G [0, T]. Recalling that, by construc- 
tion, A^- X|/ S0 <H 3 (t)-d - } = lll( x 't), and comparing (51) with the Faraday's law, we 

conclude that the electric field induced by h inside the superconductor is E(x,t) = 
w(x,t)Dp(Dh(x,t)). 

9. Further remarks 



As a consequence of the Theorem [SJ we obtain some detailed information about the 
macroscopic behavior of the penetrating magnetic field. 



Saturation time (Full penetration time). The experiments show that, if the superconductor 
is posed in an increasing external field H s , say H s (t) — > +oo as t —* +oo, after a while 
the internal magnetic field h(x,t) differs from H s (t) for a stationary amount u(x). In 
our model we obtain that if r > is such that H s (t) = supq(/io + d^), then for t > r 
h(x,t)-H s (t) = -d^(x). 

In the simpler case when H s (t) = —at, a > 0, and ho = we obtain that at time 
t = max " dn the magnetic field is penetrated in the whole superconductor and h(x,t) = 
dn(x) + H s (t) for t > r. Figures [l] and § picture the plots of the magnetic field and of the 
dissipation in the stationary configuration for anisotropic materials with a complicated 
geometry. We underline that the numerical computation based on the representation 
formula (|^) shows the presence of an hight dissipation near the points of the boundary 
with negative curvature, according to the experimental observations. 
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Figure 1 . Level sets and 3D plot of the full penetrated magnetic field and 
of the dissipation for 17 Cassini's Egg, K ellipse not centered at the origin, 
ho = 




Figure 2 . Level sets and 3D plot of the full penetrated magnetic field and 
of the dissipation for Q perturbed disk, K ellipse not centered at the origin, 
ho = 0. 

Boundary of penetration profile. If the evolution starts with ho = (no initial magnetic 
field inside f2) and the external field is nondecreasing, then at time t the induced magnetic 
field is penetrated in f2 only in the region d^{x) < H s (t)}. This gives an exact 

macroscopic description of the boundary of penetration profile, which is given by the level 
sets of the Minkowski distance from the boundary of Q. 

Hysteresis phenomenon. Assume that the evolution starts with ho = (no initial magnetic 
field inside Q) and that the external field H s (t) is increasing in [0,T], then at time T the 
induced magnetic field is h(x,T) = (—oIq(x) + H s (t)) + . Starting from this configuration, 
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Figure 3. fi perturbed disk, K ellipse centered at the origin. Profile of the 
internal magnetic field, starting from Hq = 0, during a loop of the external 
field. 



■ 




\ 









Figure 4. Magnetic hysteresis loop corresponding to the previous process. 

in [T,2T] the superconductor is subject to the external field H S (2T — t). Then we have 
h(x,2T) = min{(— d^{x) + H s (T)) + ,dn(x) + H s (0)} which in general is not zero. This 
corresponds with the experimental observation of an hysteresis phenomenon for these hard 
superconductors after an external field loop. 

On the validity of the Badid-Lopez formalism for anisotropic materials. In |Q| Fisher 
doubts the validity of the variational formulation proposed by Badfa and Lopez in the 
case of anisotropic materials. He declares that "one cannot be sure that the suggested least 
action principle chooses the Current- Voltage Characteristics correctly. Of course, if this 
principle is proved, it would be a very important tool to determine the Current-Voltage 
Characteristics via the known restriction region". Our results state that, if the power- 
law approximation for the dissipation is assumed to be reliable, then the Current -Voltage 
Characteristics induced by this approximation is the one selected by the minimization 
procedure. 
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